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Abstract 

We report a numerical study of SU(2) lattice gauge theory in the minimal Coulomb 
gauge at (3 — 2.2 and 9 different volumes. After extrapolation to infinite volume, our 
fit agrees with a lattice discretization of Gribov's formula for the transverse equal-time 
would-be physical gluon propagator, that vanishes like |k| at k = 0, whereas the free 
equal-time propagator (2|k|) _1 diverges. Our fit lends reality to a confinement scenario 
in which the would-be physical gluons leave the physical spectrum while the long-range 
Coulomb force confines color. 
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1 Introduction 



In QCD a rectangular Wilson loop W(R,T) of dimension Rx T has, asymptotically 
at large T, the form W(R,T) ~ exp[— TVw(R)], where Vw(R) is the Wilson potential. 
If dynamical quarks are present, they are polarized from the vacuum, and Vw(R) 
represents the interaction energy of a pair of mesons at separation R. In this case 
Vw(R) is not a color-confining potential, but rather a QCD analog of the van der 
Waals potential between neutral atoms. It clearly cannot serve as an order parameter 
for confinement of color in the presence of dynamical quarks, and we turn instead to 
gauge-dependent quantities to characterize color confinement. 

A particularly simple confinement scenario [|1], |2[ is available in the minimal Coulomb 
gauge.Q It attributes confinement of color to the enhancement at long range of the color- 
Coulomb potential V(R). This is the instantaneous part of the 4-4 component of the 
gluon propagator, -D 44 (x, t) = V(\x.\)5(t) + P(x., t), where the vacuum polarization term 
P(x, t) is less singular at t = 0. At the same time, the disappearance of gluons from 
the physical spectrum is manifested by the suppression at k = of the propagator 
Dij (k, fc 4 ) of 3-dimensionally transverse would-be physical gluons. This behavior is 
clearly exhibited in Fig. [l| which displays the equal-time propagators D tr (k) and D 44 (k). 

We conjecture that the color-Coulomb potential V(R) is linearly rising at large R 
(at least when asymptotic freedom holds) and that this linear rise may serve as an order 
parameter for color confinement in the presence or absence of dynamical quarks 0. 
As a practical matter, V(R) is the starting point for calculations of bound states such 
as heavy quarkonium Q. Remarkably, V(R) is a renormalization-group invariant 
in the sense that it is independent of the cut-off A and of the renormalization mass \x. 
Its Fourier transform V(k) may serve to define the running coupling constant of QCD 
by k 2 V(k) = 2 £coui(l k D- Here x o = UN/ (UN - 2JV» was calculated in §, and 
<?coui(|k|) satisfies the perturbative renormalization group equation |k|<9g cou i/<9|k| = 
/3coui(#coui)- It has the leading asymptotic behavior #c 0ul (|k|) ~ [2& ln(|k|/A coul )] _1 , 
where bo = (47r) _2 (llA — 2A/)/3, and A cou i oc Aqcd is a finite QCD mass scale. 
These formulas allow one to determine the running coupling constant of QCD from a 
numerical evaluation of the equal-time 2-point function .D 44 in the minimal Coulomb 

1 The minimal lattice Coulomb gauge is obtained by first minimizing — J2 X % r ^ r9 Ux,i with respect 
to all local gauge transformations g(x), and then minimizing — J2 X Tr 9 U x ,4 with respect to all x- 
independent but ^-dependent gauge transformations g{xi). This makes the 3-vector potential A4, 
for i = 1,2,3 transverse, diAi = 0, so A{ — . Moreover, the Coulomb gauge is the finite limit of 
renormalizable gauges ||. 
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gauge. By contrast, in Lorentz-covariant gauges, the coupling constant gv(/z) requires a 
numerical evaluation of the 3-point vertex function |J. Numerical studies of the gluon 
propagator in Landau gauge may be found in [0, |^] . 

A less intuitive but equally striking prediction concerns the 3-dimensionally trans- 
verse, would-be physical gluon propagator D i j(k., fc 4 ) = (5^ — kikj)D tT (k., fc 4 ), whose 
equal-time part is given by D tr (k) = {2tt)^ 1 J dk^D tr Ck, k^). It was proven |J as a con- 
sequence of the Gribov horizon that, for infinite spatial lattice volume L 3 , _D tr (k) van- 
ishes at k = 0, lim|k|^o -D tr (k) = 0, although the rate of approach of D tv (k, L) to 0, as a 
function of k or L, was not established, nor was it determined whether the renormalized 
gluon propagator also shares this property. This is in marked contrast to the free mass- 
less propagator (k 2 + kf)^ 1 which at equal-time is given by (27r) _1 / dk^ik 2 + fc 2 ) -1 = 
(2|k|) _1 , that diverges at k = 0. Gribov JIJ obtained, for the energy of a gluon of mo- 
mentum k, the approximate expression E(k) = |k| _1 [(k 2 ) 2 + M 4 ] 1//2 . This gives for the 
equal-time, 3-dimensionally transverse, would-be physical equal-time gluon propagator 
the approximate formula 

D tT (k) = ^(k)]- 1 = (|k|/2)[(k 2 ) 2 + M 4 ]- 1 / 2 . (1) 

In the absence of an estimate of corrections, one does not know how accurate this 
formula may be. Its accuracy, and the crucial question of the extrapolation to large L 
of D tr (\t, L) will be addressed in the fits reported here. 

We wish to confront Gribov's theory of confinement in the minimal Coulomb gauge 
[[[], [| with data from our numerical study of SU (2) lattice gauge theory, without quarks, 
in the minimal Coulomb gauge at /3 — 2.2 . Data were taken at 9 different lattice 
volumes L 4 , with L = 14, 16, 18, . . . , 30 in order to extrapolate to infinite lattice volume. 
The lattice Coulomb gauge is more easily accessible to numerical study than the Landau 
gauge because each time-slice contributes separately to the numerical average which, 
for a lattice of volume 30 4 , gives a factor of 30 gain. Details of the numerical simulation 
are described in ||12|| . (For this study we have however improved the statistics for the 
lattice volume 18 4 .) The minimal Coulomb gauge is obtained by an over- relaxation 
algorithm which, in general, leads to some one of the relative minima of the minimizing 
functional (see footnote [l|) , each of which is a different Gribov copy. As discussed in 



the Appendix of [12], the gluon propagator in Coulomb gauge is insensitive, within 



statistical uncertainty, as to which Gribov copy is attained, in agreement with similar 



2 There are corresponding results for the Landau gauge H, [|, 0. Also in Landau gauge, recent 
solutions of the Schwinger-Dyson for the gluon propagator D{k) vanish at k = [O. 
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results in the Landau gauge | 

We reported fits of V'(k) and D tT (k, L) in [12|. Our data clearly show that V(k) 
is more singular than 1/k 2 at low k, which indeed corresponds to a long-range color- 
Coulomb potential. [A linearly rising potential V(R) corresponds to V"(k) ~ l/(k 2 ) 2 ]. 
However, an extrapolation in f3 will be necessary to determine the strength of this 
singularity in the continuum limit, because U4 is quite far from the continuum limit 
at (3 = 2.2, as compared to U{, for i = 1,2, 3, due to the gauge-fixing in the minimal 
Coulomb gauge that is described in footnote |l|. In the present Letter we present a 
new fit to D tr (k, L), and most importantly, its extrapolation to infinite lattice volume 
L — > 00, using a lattice discretization of Gribov's formula. 



2 Fit to D tr (k) 

To parametrize the data we choose a fitting formula which: (i) gives a good fit to 
D tr (k, L) for all momenta k and volumes L 4 , (ii) includes Gribov's formula (j^) as a 
special case, and (iii) is physically transparent. Recall that a free continuum field with 
mass m has the equal-time propagator (27r) _1 J dk^ (A:|+k 2 +m 2 ) -1 = [2(k 2 +7n 2 ) 1 ' 2 ] -1 . 
With this in mind, consider the continuum formula Q 

D tr (k) = A~ l [z (k 2 ) a + s][(k 2 + m 2 )(k 2 + ml)}- 1 ' 2 (2) 

with singularities corresponding to poles at m\ and m 2 , which we take to be either 
a pair of real numbers or a complex-conjugate pair.[] For complex-conjugate poles 
m\ = m* 2 = x + iy, this reads 

D tT (k) = A~ l [z (k 2 ) a + s}[(k 2 + xf + y 2 )}- 1 ' 2 . (3) 

The case of a pair of real poles is obtained from this formula by taking y 2 < 0. Gribov's 
formula (ED) is recovered for a = 0.5, s = 0, x = 0, and y = M 2 (and z = 2), 
corresponding to a pair of purely imaginary poles. 

To obtain a lattice discretization of (|2|), note that for = 2sin(# M /2), the lattice 
free propagator at equal time is given by (2tt)~ 1 Jq W dO^kl+l^+ml)' 1 = (Ah\ + h\Y 1 ' 2 , 

3 In Ref. (12) we made a fit to a formula corresponding to a sum of 2 poles. However, it did not 
include the Gribov formula as a special case. 

4 According to the general principles of quantum field theory, the propagator of physical particles 
should have poles only at real positive to 2 . However, in the confined phase the gluon propagator may 
have singularities that correspond to unphysical excitations. 
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where hi = k 2 + m 2 . This suggests discretizing (g) by the substitution, 4/ii — >■ 4hi + hf, 
and similarly for 1 — > 2. For m 2 = mj 2 = x + iy, one has 4/i x + /i 2 = u + 2zyt>, where 
it = 4(k 2 + x) + (k 2 + x) 2 — y 2 , and t> = (2 + k 2 + x). This gives the lattice fitting 
formula that we shall use, 

D tr (k) = [z (k 2 ) a + s](u 2 + 42/V)" 1 / 2 , (4) 

where hi = 2sin(6 l j/2), and — 7r <&i = 2imi/L < tt, for integer rij. In the continuum 
limit one has u — > 4(k 2 + x) and v — > 2, and (^) approaches (^). 

For each lattice size L we have made a fit of the parameters z(L), s(L), ot(L), x(L) 
and y 2 {L) for the nine lattice volumes considered. Fig. Q shows our fit for D tr (h) 
for L = 18, 24, 30, and the fit to the other volumes are comparable. There is no a 
priori reason why a 2-pole fit should be accurate over the whole range of momenta 
considered.^ However, the fit is excellent for all momenta k 2 and for each L. A striking 
feature of the fit is that y 2 {L) is positive for all 9 values of L, corresponding to a pair 
of comp/ex-conjugate poles rather than a pair of real poles. Complex poles are also 



seen in the minimal Landau gauge |T3| and for gauges that interpolate between them 
15], and at finite temperature [15, 16|. By contrast, in the maximal Abelian gauge 
l|, [16|, [l7j and in Laplacian gauge poles are observed to occur at real k 2 . 



In order to do the crucial extrapolation to infinite lattice volume, we fit 5 curves^ 
z(L), s(L), a(L), x(L) and y{L) to our values of these fitting parameters. After some 
experimentation, we used the 3 fitting functions a + b/L, a + b/L + c/L 2 , and a + c/ L 2 , 
for each of the 5 parameters. Of these 3 fits, we have displayed, for each parameter, the 
one with the smallest x 2 per degree of freedom, with values reported in the captions 
of Figs 



3 Conclusions 

1. One sees clearly from Fig. ^ that the equal-time would-be physical gluon propagator 

D tT (k, L) passes through a maximum and decreases as the momentum k 2 approaches 

(for fixed L). The decrease is more pronounced as the size L of the lattice increases. 

This counter-intuitive behavior is direct evidence of the suppression of low-momentum 

components of configurations j4 tr (k) caused by the restriction to the Gribov region. 

In our simulations we consider only 3-momenta aligned along major axes &i = (0, 0, 2im/ L). Thus, 
the maximum momentum value (for each lattice side L) is 2 in lattice units, and 1.876 GeV in physical 
units (with the inverse lattice spacing a" 1 set equal to 0.938 GeV, see Ref. |l2f ). 
6 To facilitate comparison with x(L), we plot y(L) > 0. 
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2. In Fig. 3, the extrapolated value of s(oo) = 0.685 ± 0.534 is sufficiently close 
to (compared to b/L ~ 6 to 12) that we feel justified in concluding that the gluon 
propagator D tr (k) extrapolates to a value compatible with zero at k = 0. 

3. As a result, at infinite volume our fitting formula behaves like D tr (k) oc (k 2 ) a , at 
low momentum. The extrapolated value of ct(L) at infinite volume, a = 0.511 ± 0.022, 
is in striking numerical agreement with Gribov's formula ([!]), which gives a = 0.5. 

4. We have obtained an excellent 2-pole fit for D tr (k, L). Our fit indicates that the 
poles occur at complex m 2 = x(L) ± iy(L). The real part x(L) extrapolates to 0.006 ± 
0.024 at L = oo. Remarkably, a pair of poles in k 2 at purely imaginary m 2 = ± iy 
agrees with the Gribov equal-time propagator D tr (k) = (2|k|) _1 [l + M 4 (k 2 ) -2 ] -1 / 2 , 
with y 2 = M 4 . Moreover at large k, this gives a leading correction to the free equal- 
time propagator of relative order (k 2 )~ 2 with coefficient of dimension (mass) 4 . It may 
not be a coincidence that this is the dimension of the gluon condensate (F 2 ), which is 
the lowest dimensional condensate in QCD. Because of the theoretical suggestiveness 
of our result, we are encouraged to report the values m 2 = ± iy, for y = 0.375 ±0. 162 
in lattice units, or y = 0.330 ± 0.142 GeV 2 , M = y 1 ' 2 = 0.575 ± 0.124 GeV for the 
location of the gluon poles in k 2 . Here it is assumed that we are already in the scaling 
region for D tr (k), which remains to be verified by further numerical studies. That this 
is not entirely unreasonable is suggested by the fact that scaling for the Landau gauge 
propagator in SU(2) has been observed || in the range 2.1 < /3 < 2.6. Note also that 
Gribov derived his formula in the continuum theory, namely at infinite (3, so on the 
theoretical side it is to be expected that our fit will remain valid throughout the scaling 
region. 

5. The observed strong enhancement fl2| of the instantaneous color-Coulomb po- 
tential V(k) and the observed suppression of the equal-time would-be physical gluon 
propagator D tr (k) both at low k, strongly support the confinement scenario of Gribov 
[|TJ In addition to this qualitative agreement, we note excellent numerical agreement 
of our fit to the formula D tT (k) = ,z(|k|/2)[(k 2 ) 2 + M 4 ]" 1 / 2 . Although this formula can- 
not be exact, it appears that deviations from it are relatively weak at both high- and 
low-momentum regimes. If this excellent fit is maintained at larger f3 values, then it 
appears that we have obtained a quantitative understanding of D tr (k). 
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Figure 1: Plot of the gluon propagators D tr (k) and .D 44 (k) as a function of the square 
of the lattice momentum k 2 for L = 28 (symbols * and O respectively) and L = 30 
(symbols A and y respectively). Notice the logarithmic scale in the ?/ axis. Error bars 
are one standard deviation. 
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Figure 2: Plot of the gluon propagator D tr (k) as a function of the square of the lattice 
momentum k 2 for L = 18(*), 24(A) and 30(O) an d fits of these data using Eq. (^). 
Error bars are one standard deviation. 
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Figure 3: Fits of the parameters z(L)(*) and s(L)(A) plotted against 1/L. For 
z(L) =a + b/L, we obtain a = 7 A3 ± 0.36, b = -72.5 ± 7.7, and xVd.o.f = 1.13. For 
s(L) =a + b/L, we obtain a = 0.685 ± 0.534, b = 165.6 ± 12.7, and x 2 /d.o.f = 0.183. 
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Figure 4: Fits of the parameters x(L)(*) and y(L)(A) plotted against 1/L. For 
x(L) = a + c/L 2 we obtain a = 0.006 ±0.024, c = 59.0 ±9.9, and x7 d -°- f - = 0.858. For 
y(L) = a + b/L + c/L 2 we obtain a = 0.375 ± 0.162, b = 13.7 ± 6.3, c = -101.6 ± 61.0, 
and x7 d -°- f - = 1-024. 
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